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Abstract 

We update our estimates of charged and neutral current neutrino total cross sections on 
isoscalar nucleons at ultrahigh energies using a global (x, Q 2 ) fit, motivated by the Froissart 
bound 1 ] , to the F2 (electron- proton) structure function utilizing the most recent analysis [2] of 
the complete ZEUS and HI data sets from HERA I. Using the large Q 2 , small Bjorken-x limits 
of the "wee" parton model, we connect the ultrahigh energy neutrino cross sections directly to 
the large Q 2 , small x extrapolation of our new fit, which we assume saturates the Froissart 
bound pp. We compare both to our previous work [3J, which utilized only the smaller ZEUS 
[I] data set, as well as to recent results \5\ of a calculation using the ZEUS-S based global 
perturbative QCD parton distributions using the combined HERA I results as input. Our new 
results substantiate our previous conclusions 3 , again predicting significantly smaller cross 
sections than those predicted by extrapolating pQCD calculations to neutrino energies above 
10 9 GeV. 

1 Introduction 

Large experiments seeking evidence of ultrahigh energy (UHE) cosmic neutrinos have been going 
on in earnest for more than a decade^] Searches have been underway with detectors scanning for 
UHE cosmic neutrino induced events in large volumes of water [6j [7] , ice [8j [9j [TUJ [11] , Earth's 
atmosphere [TH and the lunar regolith [TJ]. No clear indication of a cosmic neutrino event 
has yet been reported, but by assuming extrapolations of neutrino-nucleon cross sections from low 
energy data to UHE, the experiments have all reported bounds on UHE neutrino flux models. Put 
together, the experimental bounds on neutrino fluxes now cover energies up to 10 GeV. 

Given the importance of neutrino cross section estimates in the process of interpreting data 
(i.e., placing flux bounds) or planning new experiments (i.e., EUSO [IS])) having a range of realistic 
predictions of neutrino cross sections available is crucial to the ultimate success of the UHE neutrino 



1 For purposes of this paper, we take UHE to mean E v > 10 6 GeV and cosmic to mean origins that are galactic, 
extragalactic or cosmogcnic. 
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physics enterprise. It was in this spirit we applied an extrapolation to UHE of a global fit to F2 
structure function data [16] from ZEUS [4] to an estimate of the UHE neutrino-nucleon total cross 
section [3]. The fit function was well motivated, invoking the Froissart hi 2 (s) asymptotic bound on 
hadronic cross sections pQ, |17j . and with only a handful of parameters, it provided an excellent fit 
to the ZEUS data set which covered an impressive range of Q 2 and x values. Though the agreement 
with the conventional perturbative quantum chromodynamics (pQCD) prediction was good in the 
range of E v where the cross section is dominated by the HERA (x, Q 2 ) domain, when extrapolated 
to energies above 10 9 GeV our prediction fell well below the pQCD extrapolation, having important 
implications for the estimate of event rates for the highest energy experiments. 

With the significance of these issues in mind, we devote this paper to updating the results of 
our earlier calculation by using the recently published joint ZEUS and HI analysis of all the HERA 
I data [2]. In particular, this new analysis reconciles some tension between relative normalization 
between the two experiments, making a complete evaluation of errors for the combined data set. 
We fit the new results as before, recalculate <j„n{E v ) and compare both with our previous results 
[3], as well as new results of a standard pQCD calculation based on a recent determination of quark 
distribution functions [5] using a earlier analysis [13], updated to include all the HERA I data. We 
confirm our previous conclusion that our unitarity bounded cross sections predict values that are 
significantly lower than those found from unscreened pQCD at neutrino energies above 10 9 GeV. 
Moreover, even when we extrapolate our fits many orders of magnitude in energy above the region 
covered by the HERA data, the greatly expanded and re-analyzed data base leads to only very 
modest changes in our new values for <j v n when compared to our previous ones [3]. 



2 The UHE v — N cross section from the full HERA I data 
set 

We assume two full families of active quarks. In the "wee parton" picture, assuming that valence 
quark contributions are small at UHE and that sea quarks all contribute the same (equipartition 
of flavors), the neutrino-isoscalar nucleon N (N=(p+n)/2) double differential cross section can be 
written as 

d 2 a v _ 2G%mE v ( M 2 



dxdy 



{Ev) = 2G^ (q^t) >< Mi(x, Q 2 )] (1 + (1 - y)% (1) 

where E v is the neutrino energy in the rest frame of the nucleon. The electron-proton structure 
function F 2 (x, Q 2 ) is provided by our fit to the combined ZEUS and HI analysis of the full HERA 
I data set [2J, described in the following section. The label V takes the values V = CC (charged 
current) and NC (neutral current), where Mcc=M\y, M^c—Mz, m is the mass of the nucleon, Gf 
is the Fermi weak coupling constant, and kcc—jq and ktvc=0-298 for two full families, using the 
2008 PDG [19] value of sm 2 (9 w ) at Q 2 =M| in the WS scheme. The factors relating the electron- 
proton Fi structure function (found in deep inelastic scattering) to the CC and NC equivalents, 
using the '"wee parton" model, are derived in [3] . The general quark parton CC and NC expressions 
are given in [20] . 
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2.1 Kinematics 



In Eq.([T]), the inclusive double differential cross section for the CC (NC) neutrino inclusive cross 
sections for the processes ve(k)+N(p) — > £(k')+Xcc [ve(k')+Xjfc)) ^ = e, /i, r, depends upon the 
invariants s — (k+p) 2 , —Q 2 = (k — k') 2 , and p- (k — k 1 ). The scaling variables are x — Q 2 /2p- (k — k 1 ) 
and y — p-(k — k')/p-k. If E v 3> m, s — 2mE v and Q 2 = 2mE v xy. In the rest frame of the nucleon, 
y = 1 — E' u /E u , the fractional energy transferred from the incoming neutrino to the outgoing lepton 
or neutrino. 

In the integrations over x and y, whose results we will report below, we set a limit to the 
minimum value of Q 2 , which avoids possible numerical problems with singular behavior of the 
integrands as x or y goes to zero. Given the relationship shown in the preceding paragraph, we see 
that this amounts to limiting x and y to the ranges Q 2 nin /(2mE v ) < x < 1 and, taking y as the 
initial integration, Q 2 nin j '{2mE l/ x) < y < 1. 

As recognized early in |5I] . the vector boson propagator factor (M V /(M 2 + Q 2 )) 2 acts to cut 
off the integrand for Q 2 > My, effectively selecting a range of small x — reaching somewhat below 
x ~ M v / (2mE u ) — that makes substantial contributions to the total cross section. For the range of 
neutrino energies we consider in this work, 10 4 GeV < E v < 10 14 GeV, this means that we probe 
x values in the range 0.1 < x < 10 -11 ! 

2.2 Global fit to combined ZEUS and HI F 2 (x,Q 2 ) results 

We follow the work of Berger, Block and Tan [16], who obtained a remarkably good fit to ZEUS 
small x results for F 2 (x, Q 2 ) [4], using a 6 parameter model guided by empirical analysis [17] 
of experimental data on high energy hadronic cross sections that empirically demonstrated the 
saturation of the In 2 s Froissart bound [I]. They assumed that Deep Inelastic Scattering (DIS) 
was 7*p hadronic scattering and applied the In 2 s Froissart bound saturation successfully to the 
F2(x, Q 2 ) data [3]. For Q 2 m 2 , this bound in s translates into a ln 2 (l/a;) bound on the small x 
behavior of F 2 {x, Q 2 ). We here apply the same fit function to the recently reported results from a 
combined ZEUS and HI determination of DIS e^p cross sections [2] , which provides very accurate 
values of F2(x,Q 2 ) over a large region of the (x, Q 2 ) plane. Using Eq. (JT]), we then predict UHE 
neutrino cross sections based on the extrapolation of the analytic expression F 2 into the Q 2 > M§ 
and x < 1CU 10 regions needed to evaluate the cross sections up to 10 12 GeV and beyond. 
The global fit function takes the form, for small x, 

F 2 {x,Q 2 ) = (l-xX-^+AiQ^H—^—^-] 

1 — Xp XL — Xp 

+B{Q 2 )\A—±—^]), (2) 

X 1 — Xp 

where 

A(Q 2 ) = a Q + ai In Q 2 + a 2 In 2 Q 2 

B(Q 2 ) = b + b 1 \nQ 2 + b 2 \n 2 Q 2 . (3) 

A x 2 minimization global fit of Eq.@ and Eq.Q was made to all of the HERA combined experi- 
mental F 2 data in both Q 2 and a;, with x < xp = 0.11. At the point x P = 0.11, dF 2 (xp, Q 2 )/dQ 2 =0 
for all Q 2 . This procedure differs slightly from that in [IB], where Fp was fixed at the value 0.41; 
as in [16| . we again used the "Sieve" algorithm [55] to eliminate "outlier" datum points, using a 
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^Xmax = 6 cut. Fp, the value of F 2 at xp, along with the other 6 parameters, together with their 
errors, are listed in Table 1. Also shown are the renormalized minimized x 2 value |22j . the number 
of degrees of freedom and the renormalized \ 2 P er degree of freedom for our new analytic form for 
the combined ZEUS and HI results [2]. 

The large x ~ 1 region, which contributes very little the UHE cross sections, is fitted by the 
form 

Table 1: Results of a 7-parameter fit to the HERA combined data for F 2 (x,Q 2 ) for 0.85 < Q 2 < 3000 
GeV 2 and x < 0.1. The 

Xmin is renormalized by the factor 1Z to take into account the effects of the cut at 



Parameters 


Values 


a 


-8.471 x 10" 2 ±2.62 x 10" 


-3 


ai 


4.190 x 10~ 2 ± 1.56 x 10" 


-3 


a 2 


-3.976 x 10- 3 ±2.13x 10" 


-4 


bo 


1.292 x 10" 2 ±3.62 x 10" 


-4 


h 


2.473 x 10" 4 ± 2.46 x 10 


-4 


b 2 


1.642 x 10- 3 ±5.52x 10 


-5 


F P 


0.413 ±0.003 


Y 2 ■ 


352.751 




^ X Xmin 


391.377 




di. 


335 




^ X Xmin/ di - 


1.17 





Equations © and ((H) obviously guarantee continuity of F2(x,Q 2 ) at x = xp for all x. Requiring 
that dF 2 (x,Q 2 )/dx is also continuous at x — xp determines p{Q 2 ) in terms of xp 7 Fp and A(Q 2 ), 
Eq. ©. 

3 Total charged and neutral current neutrino-nucleon cross 
sections 

The total deep inelastic scattering neutrino-nucleon cross section is obtained by the double integral 
of Eq.(fT|) over Q 2 min j (x2mE v ) < y < 1 and Q 2 min l '{2mE y ) < x < 1. Tables 2 and 3, calculated using 
Qmin = 0-01 GeV 2 , show the charged and neutral current total i/-isoscalar nucleon cross section 
values in cm 2 , at energies from 10 4 to 10 14 GeV. The column labels cthera, o'bbmt, cc-ss and 
&GQRS refer, respectively, to a) the calculations based on our new fit-function — Eqs. ([2]) and (|4]) — 
to the combined ZEUS and HI measurements of F 2 [2], b) the published Berger, Block, McKay 
and Tan (BBMT) values based on the same fit function applied to the older ZEUS F 2 values [3], c) 
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a recent calculation [5] using ZEUS-S [TB] parton distribution function (pdf) fits with NLO pQCD 
evolution [23], and d) the results from the 1998 paper of Gandhi, Quigg, Reno and Sarcevic [20] 
using CTEQ4 fits with DIS pQCD evolution. 

Table 2: Charged current cross sections, in cm 2 , as a function of E v , the laboratory energy in GeV; 
&hera, from our model using a fit to new HERA results; a b bmt, from our earlier fit to ZEUS 
data [3j <Jc—SSi from a standard evaluation from recent pdfs [5] ; ogqrs, from the 1998, CTEQ4 
based pdf [2D] , Note that the latter two only evaluate values up to E v = 10 12 GeV. 



E v (GeV) 


&hera( 


cm 2 ) 


<JBBMT(cm 2 


) 


<?c-ss (cm 2 ) 


&GQRs{ 


cm 2 ) 


10 4 


4.57 10- 


35 


3.80 10" ab 




4.50 10" 35 


4.62 10- 


-35 


10 5 


2.11 10- 


34 


1.91 10~ 34 




1.95 10~ 34 


2.02 10- 


-34 


10 6 


6.88 10- 


34 


6.87 10- 34 




6.01 10~ 34 


6.34 10- 


-34 


10 7 


1.90 10- 


33 


1.94 10~ 33 




1.60 10~ 33 


1.75 10- 


-33 


10 s 


4.48 10- 


33 


4.49 10~ 33 




3.87 10~ 33 


4.44 10- 


-33 


10 9 


9.09 10- 


33 


8.90 10" 33 




8.76 10- 33 


1.05 10- 


-32 


10 10 


1.64 10- 


32 


1.58 10~ 32 




1.87 10~ 32 


2.38 10- 


-32 


10 11 


2.72 10- 


32 


2.57 10" 32 




3.81 10~ 32 


5.36 10- 


-32 


10 12 


4.21 10- 


32 


3.92 10~ 32 




7.40 10~ 32 


1.18 10- 


-31 


10 13 


6.17 10- 


32 


5.68 10~ 32 










10 14 


8.68 10- 


32 


7.92 10~ 32 











The results are summarized graphically for the evaluations (Jhera and ac-SS for CC, together 
with cthera for NC in Fig. [1] Our new results are shown as a (red) solid upper curve for ace 
and a (blue) dot-dashed lower curve for <jnc- The perturbative QCD results of [5] for acc are 
shown by the (black) dashed curve up to E v = 10 12 GeV, the maximum energy they report. The 
latter reference does not present NC results. Though they are not shown in Tables 2 and 3, we 
note in passing that the errors in our cross section values (including correlation terms) that are 
generated by the uncertainties in the fitted parameters range between 1% and 1.7 % in the interval 
from 10 5 GeV to 10 14 GeV, and are thus quite small in comparison to the differences between our 
calculations and the pQCD calculations for energies above 10 9 GeV, as illustrated in Fig. 1. In both 
the CC and the NC cases, we see from the tables and from the plot that the pQCD extrapolation 
for g v n{E) rises significantly above our "Froissart" extrapolation for E v > 10 9 GeV. 

4 Summary and Conclusions 

We have updated our previous prediction of the ultrahigh energy CC and NC neutrino cross sections 
by repeating our earlier calculation with a recently published, combined ZEUS and HI data set 
with nearly double the statistics compared to that used in our calculation based on the 2001 ZEUS 
results alone. Besides increasing the size of the data sample, the new analysis made a thorough 
study of the correlated errors, resolving the tension between the individual analyses. The results 
we report here use a 7 parameter fit, achieving a x 2 /d.f. of 1.17, similar to that obtained by the 
previous 6 parameter fit to the much smaller ZEUS data set. As indicated by the tabular and 
graphical results summarized in Sec. 3, the CC cross section predictions of pQCD and our model 
in the UHE range 10 6 <E V < 10 8 GeV agree reasonably well. This is not too surprising, since the 
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Table 3: Neutral current cross sections, in cm 2 , as a function of E„, the laboratory energy in GcV; 
o~hera, from our model using a fit to new HERA results; o~bbmt, from our earlier fit to ZEUS 
data, [3J; (Jgqrs, from the 1998 CTEQ4 based pQCD extrapolation of Ref. {2"U1 . which quotes 
values only up to 10 12 GeV. Note that Ref. [5] does not report NC results. 



E v (GeV) 




obbmt (cm 2 ) 


&GQRS (cm 2 ) 


10 4 


1.60 10" 35 


1.32 10~ 35 


1.58 10~ 35 


10 5 


7.88 10~ 35 


7.03 10" 35 


7.67 10" 35 


10 6 


2.65 10~ 34 


2.65 10~ 34 


2.60 10" 34 


10 7 


7.55 10" 34 


7.74 10" 34 


7.48 10" 34 


10 8 


1.82 10~ 33 


1.83 10~ 33 


1.94 10" 33 


10 9 


3.76 10" 33 


3.70 10~ 33 


4.64 10" 33 


10 10 


6.89 10~ 33 


6.63 10" 33 


1.07 10" 32 


10 11 


1.15 10" 32 


1.09 10~ 32 


2.38 10~ 32 


10 12 


1.79 10~ 32 


1.67 10~ 32 


5.20 10~ 32 


10 13 


2.64 10" 32 


2.44 10~ 32 




10 14 


3.73 10~ 32 


3.40 10~ 32 





average x and Q 2 values contributing to a v (E ll ) in this range are both consistent with the "wee 
parton" model we invoke and lie within the HERA data range. What is somewhat surprising is the 
good agreement for E v <10 6 GeV. 

When E v > 10 9 GeV, the ln 2 (l/ic) term dominates the values of F%, and the pQCD-based value 
of (Jvn{E u ), which is growing approximately as a fractional power of x, become clearly larger than 
our fit based on "Froissart-unitarity" growth with E v . Note, however, that [5] finds that HERA I 
data — the input to the ZEUS-S global pdf fits [18] — indicate that the growth of the cross section is 
not quite a power law and falls a bit below the earlier determination by Gandhi, Quigg, Reno and 
Sarcevic (GQRS) at energies above 10 9 GeV. At the value E v = 10 12 GeV, for example, the 
Ref. 20 value for the CC cross section is a factor of 1.5 larger than that of Ref. [5]. 

Though the data set used for this work is much larger, basically doubling the statistics, and the 
error analysis more complete than the set used in our previous study [3], there is qualitative and 
general quantitative agreement between the two at the very highest energies, even when extrapolated 
far above the kinematic region where the data lie. We view this as evidence for the stability and 
robustness of our extension of the data to ultrahigh energies. 

In conclusion, our results based on the complete HERA I data set confirm the message of our 
previous work that hints of unitarity constraints may already be present in the HERA deep inelastic 
scattering evaluations of Fi. The unitarity constraints are most evident in the predictions of the 
neutrino-nucleon total cross sections in the experimentally important realm of the highest energy 
cosmic rays, E v >10 9 GeV, where cosmic ray experiments [H] [T3] and others such as RICE [5], 
ANITA [TT]. GLUE [14] and FORTE [10] have set limits on the fluxes of cosmic neutrinos. The 
cross-section issue will ultimately be settled by an analysis such as that outlined in Ref. [24] . 
There the experimental requirements are laid out that would enable one to distinguish among the 
unscreened QCD (the pQCD approach), screened QCD, estimated in the "dipole" approximation 
(whose predictions resemble ours at the highest energies), and a radical non-linear picture like the 
quark-gluon condensate. Meanwhile, in planning for new experiments and analysis of ongoing ones, 
our direct extrapolation from data, which does not use QCD evolution [23] but builds in unitarity 
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10 -35 Ll 

10 4 10 6 10 8 10 10 10 12 10 14 

E(GeV) 

Figure 1: This plot shows the total charged current and neutral current cross sections for neutrino 
scattering on an isoscalar nucleon target, with E v in GeV and a v in cm 2 . For charged currents, 
HERA, our fit, is the (red) solid curve and the next-to-leading order C-SS [5] fit is the (black) 
dashed curve. For neutral currents, the (blue) lower dot-dashed curve is for our new HERA results. 
Note that neutral currents are not included in the C-SS study. 

constraints, should play a critical role in estimating event rates and detection thresholds. 
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